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The conventional duality analysis is employed to identify a location of a critical point on a uniform 
lattice without any disorder in its structure. In the present study, we deal with the random planar 
lattice, which consists of the randomized structure based on the square lattice. We introduce 
the uniformly random modification by the bond dilution and contraction on a part of the unit 
square. The random planar lattice includes the triangular and hexagonal lattices in extreme cases 
of a parameter to control the structure. The duality analysis in a modern fashion with real-space 
renormalization is found to be available for estimating the location of the critical points with wide 
range of the randomness parameter. As a simple testbed, we demonstrate that our method indeed 
gives several critical points for the cases of the Ising and Potts models, and the bond-percolation 
thresholds on the random planar lattice. Our method leads to not only such an extension of the 
duality analyses on the classical statistical mechanics but also a fascinating result associated with 
optimal error thresholds for a class of quantum error correction code, the surface code on the random 
planar lattice, which known as a skillful technique to protect the quantum state. 



I. INTRODUCTION 

The duality, which bridges two distinct concepts or 
situations, is often found and plays important roles in 
various realms of physics and mathematics. We focus on 
the duality in statistical mechanics, which is embedded 
in the symmetry of the partition function at low and high 
temperatures. This symmetry allows us to identify the 
locations of the critical points for various spin models 
such as the Ising and Potts models 0, [1] . 

The most remarkable development on the duality in 
the spin systems is the successful application to the ran- 
dom spin systems, such as the spin glass, which has phase 
diagrams of complicated structures and exhibits rich be- 
haviors. Since the series of work revealed that the 
duality is capable to estimate the approximate values of 
the critical points even for such complicated systems, the 
duality has been developed as one of the rare success- 
ful approaches to finite-dimensional spin glasses, where 
very little analytical work exists. More precise analysis 
of the locations of the critical points in spin glasses has 
been performed by employing a real-space renormaliza- 
tion technique, i.e., the partial summation of the degree 
of freedom 

Very recently, the duality also plays an important role 
to provide a powerful tool in other fields. For instance, 
it has opened a way to analyze the theoretical limita- 
tion in the classical/quantum coding theory, where order- 
disorder phase transition corresponds to the theoretical 
limitation of the error correction codes. Specifically, the 
optimal error thresholds of a class of quantum error cor- 
rection codes, so-called surface codes [9(, have been eval- 
uated in terms of the duality with the real-space renor- 
malization [l(J HH, On the other hand, newly developed 
quantum error correction codes (l2j also motivate us to 
study the associated spin glass models, which have not 
been considered so far. 



In the present paper, we extend the conventional dual- 
ity analysis to a broader class of the lattices, which have 
various types of randomness. Specifically, we deal with 
random planar lattices generated by random modifica- 
tion, bond dilution and contraction. The random planar 
lattice consists of a mixture of the triangular and the 
hexagonal lattices through a single parameter to control 
its structure. Since the conventional duality analysis is 
limited to a few self-dual lattices (e.g. the square lattice), 
we cannot deal with such a random planar lattice on all 
parameter values systematically. In order to overcome 
this situation, we establish a unified and systematic ap- 
proach to analyze the random planar lattices by use of 
the duality with the real-space renormalization. In cer- 
tain extreme cases, we can reproduce the critical points 
on the triangular and hexagonal lattices, where the so- 
called star-triangle transformation works well with the 
duality analysis [13|,[l4]j. In this sense, the present frame- 
work can be regarded as a natural extension of the star- 
triangle transformation to a broader range of the random 
planar lattices. 

II. DUALITY ANALYSIS 

The structure of the random planar lattice is based on 
the square lattice, where we can obtain the exact location 
of the critical points via the duality analysis by virtue 
of its special symmetry called self-duality [l], [|[ . Let us 
here briefly review the duality analysis by taking a simple 
case. We consider the q-state Potts spin system with the 
following Hamiltonian 

H = -^2j ij S{S it S j ), (1) 

m 

where S(x,y) is Kronecker's delta and Si takes 
Q, ,q — 1. Now, we set the couplings to 
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be a positive value J > for simplicity. The du- 
ality is a symmetry argument in the low and high- 
temperature expansions of the partition function Z = 
Hsi Yluj) exp((3J5(Si, Sj)) [1]. From a modern point of 
view, these expansions can be related with each other 
by the g-component discrete Fourier transformation of 
the local part of the Boltzmann factor, namely edge 
Boltzmann factor Xk — exp(/3 J6(k, 0)) for a single bond 
0. Each term in the low-temperature expansion can 
be expressed by x^. On the other hand, the high- 
temperature expansion is expressed by the dual edge 
Boltzmann factor through the discrete Fourier transfor- 
mation x\ = J2t x k exp(i2ivkl/q)/y/q. By using this fact 
we can obtain a double expression of the partition func- 
tion as 

(xo) Nb z(u u U2, ■■■) = (x*) Nb z{uI, «£,.-. ), (2) 

where z is the normalized partition function 
z(ui,u 2 , ■ ■ ■ ) = Z/(x ) Nb and z(u\, m|) " * ' ) = 
Z/(xq) Nb . We here define the relative Boltzmann 
factors Uk = x^/xq = exp(— j3J) and u* k — x* k /x$ — 
(exp(/3J) - l)/(exp(/3J) + q- 1). The well-known duality 
relation exp(/3*J) = (exp(/3J) + q- l)/(exp(/3J) - 1) 
is obtained by rewriting u k by Uk with the different 
coupling constant, which implies the transforma- 
tion of the temperature. Its fixed point condition 
exp(^ c J) = (exp(/3 c J) + q — l)/(exp(/3 c J) - 1) gives 
the critical temperature. We can evaluate the exact 
value of the internal energy and the rigorous bound for 
the specific heat at the critical point from the relation 
of the partition function ^ through the appropriate 
derivatives as well as its location (l5j . 

If we are interested in derivation of the location of 
the critical point, we can use an alternative way. We 
impose that both of the principal Boltzmann factors Xq 
and Xq with edge spins parallel are equal, (i.e. xq = Xq). 
We reproduce the equation exp(/? c J) = (exp(/3 c J) + q — 
l)/(exp(/? c J) — 1). Notice that the symmetry of the low- 
and high-temperature expansions on the partition func- 
tions is embedded compactly in the edge Boltzmann fac- 
tors only on the single bond. In particular, the above 
argument is extremely simplified if we use the principal 
Boltzmann factors. 

In the duality analysis, two expansions are related in- 
between the original and its dual lattices. In the dual 
lattice, the site and plaquette of the original lattice are 
exchanged with each other, and the bonds are perpen- 
dicular to those on the original lattice. The bonds de- 
scribed on the dual lattice are perpendicular to the orig- 
inal bonds. The square lattice is one of the self-dual 
lattices, which share the same structure even after the 
transformation. Usually, the duality analyses are limited 
to the case on the self-dual lattice. For non-self-dual lat- 
tices, such as the triangular and hexagonal lattices, the 
duality analysis gives only the relationship between crit- 
ical points for mutually dual pair 0, EH • 

In the case we will consider, the coupling constants Jy 
are multivariated due to the existence of the quenched 



randomness on the interaction and/or structure of the 
lattice, which we will introduce later. This hampers the 
straightforward application of the conventional duality 
analysis. We then employ the replica method, which gen- 
eralizes the duality analysis to deal with the quenched 
randomness [1,0]. 

Let us consider the duality for the replicated partition 
function as [Z n ] (written simply as Z n below), where 
[• • • ] is the configurational average for the quenched ran- 
domness according to the distribution functions of J+j. 
The (/"-component discrete Fourier transformation again 
leads us to the double expressions of the replicated par- 
tition function as 

Z n (x ,xi,- ■ ■) = Z n (x* ,xl, ■ ■ ■ ). (3) 

We extract the principal Boltzmann factors from both 
side of Eq. ([3]), and obtain 

(x Q ) NB z n (u 1 ,u 2 , ■■■) = (xq) Nb z n (ui,u 2 , ■■■), (4) 

where Z n (ui, u 2 , ■ ■ • ) = Z n (xo) NB and z n (u*i, uJj, • • • ) = 
Z n ]/ {xq) Nb . In the case with the quenched randomness, 
unfortunately we cannot replace u* k by Ufc as the pre- 
vious case without the quenched randomness, since the 
replicated partition function is multivariable in general. 
Nevertheless, we assume that Xq is equal to Xq again as 
done for the case without the quenched randomness, and 
take the limit of n — > following the recipe of the replica 
method. As a consequence, we estimate the approximate 
location of the critical points for spin glasses. 

By taking a wider range of the local part of the Boltz- 
mann factor rather than a single bond, we can improve 
the estimates, and the resultant values are expected to 
asymptotically reach the exact solutions 0, @- For in- 
stance, in the case on the square lattice, we define the 

is) 

cluster Boltzmann factor x k , where the subscript k de- 
notes the configuration of the edge (white-colored) spins, 
and s distinguishes the size of the cluster as in Fig. [T] 
Then we can rewrite the equality from the conventional 
duality ((3]) as 

Z n (xQ \x\ , ■ ■ ■ ) = Z{xq , Xi , • • • ). (5) 

Again we extract the renormalized-principal Boltzmann 
factors x^ and x^^ from both sides. 

= (x^Y B, 4< (s \u 2 is \---), (6) 

where z^ is the normalized partition function de- 
fined as z^(u[ s \u { 2 s \---) = Z n /{x[ s) ) N B\ and 
z{ul {s \u 2 (s \---) = Z n /(xo {s) ) N B\ We here define the 
renormalized-relative Boltzmann factors — x^ /x^ 

and u k ^ — x*£ s ' /x*q S \ We set the equation to esti- 
mate the more accurate location of the critical point than 
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FIG. 1: The unit cluster for the square lattice with four bonds 
s = 1 (left upper panel). The left lower panel describes the 
single bond once we considered for the simple application of 
the conventional duality with the replica method. The right 
panel denotes the larger cluster with 16 bonds s = 2 to iden- 
tify more precise location of the critical points. The black 
circles denote the spin variables summed over. The white 
circles express the spin variables without summation. 



as inspired by the previous case without the 



Xq = 

quenched randomness, 



» 



(7) 



where Xq and XQ (b> express the renormalized Boltzmann 
factor with the edge spins parallel as the detailed calcu- 
lation shown below. The conventional duality with the 
single bond can be understood as the special case s = 0. 

Although the above method is not exact, it is believed 
that the estimations systematically approach the exact 
solutions for the critical points of spin glasses by increas- 
ing the size of the cluster 0, H, [JJ[ . 
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FIG. 2: Structure of the random planar lattice. The lower 
panels depict a part of the unit cell of the whole lattice. The 
two types of modification, "bond dilution" and "bond con- 
traction" are assigned on each dashed bond. The right panels 
give the dual lattice, while the left ones express the origi- 
nal lattice. Readers can confirm the self-duality holds on the 
structure of the lattice even with existence of a part of the 
modifications. 




III. RANDOM PLANAR LATTICE 

Let us take a step toward the main issue in the present 
study. We introduce two types of quenched randomness 
for the coupling constant Jij to form the random planar 
lattice. One is the bond dilution, which means that we 
remove the bond. The other is bond contraction by which 
we put two of the sites on the end of the bond together. 
The bond dilution and contraction expressed by Jy = 
and Jij — oo are given with the probability p m ix and 
1 — Pmix, respectively. We chose a single bond on each unit 
cell with such randomness as in Fig. [5] We describe the 
representative instance of the random planar lattice with 
Pmix = 1/2 in Fig. [3] If we take the extreme limit p m i x — > 
0, the resultant lattice becomes the triangular lattice as 
in Fig. 0] (left). On the other hand, when p m i x — > 1, it 
becomes the hexagonal lattice as in Fig. [5] (left). These 
two lattices are known as a mutually dual pair. While, 
in this case, we can obtain a relation between two critical 
points on the mutually dual pair of lattices, we cannot 
identify the location of the critical point. Thus we needed 
to employ an additional tool to derive the location of the 
critical point, such as the star-triangle transformation 
[l3l Il4j |. However, once we regard the random planar 



FIG. 3: One of the realizations of random planar lattice (left), 
and its dual one (right). 




FIG. 4: Extreme limit of p m ix. = 0. The original lattice in the 
left panel becomes the triangular lattice. On the other hand, 
the dual one in the right-hand side is the hexagonal lattice. 
For guide to eyes, we rearrange the location of the sites as the 
modifications of the unit cell on the lower panels. 
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FIG. 5: Extreme limit of p m i K . = 1. The original lattice in 
the left panels turns into the hexagonal lattice, while the dual 
one in the right-hand side becomes the triangular lattice. The 
same symbols are used as those in Fig. [4] 



lattice as the square lattice with quenched randomness 
such as Jij = and Jij = oo given by probability p m i x 
and 1 — Pmix, respectively, we may perform the duality 
analysis on such a random spin model by virtue of the 
duality of the original square lattice. Actually it will be 
shown below that the self-duality on the original square 
lattice always holds for any p m i X . This is the reason why 
we can employ the duality analysis with the real-space 
renormalization on the cluster of the unit cell as in Fig. Q] 
to analyze the location of the critical point on the random 
planar lattice. 

Beyond the present case, we can generalize the condi- 
tions on the class of random planar lattices, where the 
duality works well by considering modification starting 
from the square lattice. The key point is to hold the self- 
duality of the square lattice even for the location of the 
quenched randomness, i.e., bond dilution and contrac- 
tion. As in Fig. [2 the dual lattice has the same struc- 
ture including the locations of the bonds to be modified 
as that of the original one. One can confirm this fact 
by rotating the dual lattice by 90 degrees. Following 
this rule, we demonstrate other random planar lattices 
by tiling the unit cells as in Fig. [B] The first instance 
as in Fig. [5] (left) is trivial as a spin model, since it does 
not undergo any phase transition in finite temperature. 
The other as in Fig. [S] (right) possesses a complicated 
structure, on which the spin models would show finite- 
temperature phase transitions. 

One might imagine the generalization of the random 
planar lattice enables us to investigate most of the planar 
lattices by considering appropriate modifications. How- 
ever, in most cases, one cannot find such appropriate 
modifications. The system on the random planar lat- 
tice possesses inhomogenous interactions, if we regard the 
bond under dilution or contraction as the different inter- 
action denoted by Jij — or oo from the uniform ones 
Jij = J. Let us take a counterexample here. We consider 
a random planar lattice by use of the unit cell as shown 



FIG. 6: Other random planar lattices where one can perform 
the duality analysis. The left panel has a simple case by tiling 
of the unit cell with four bonds. The right panel is a more 
complicated case by the unit cell with 16 bonds. 




FIG. 7: A random planar lattice where one can not obtain the 
critical points. In the lower left panel, we describe the unit 
cell. The extreme limit p m i x = gives the (4, 8 2 ) lattices. This 
is the notation of the uniform tilings, in which all polygons are 
regular and each vertex is surrounded by the same sequence 
of polygons. The notation (4, 8 2 ) expresses that every vertex 
is surrounded by one square and two octagons. 



in Fig. [7l The extreme limit p m j x = becomes the (4, 8 2 ) 
lattice. In contrast to the successful instance, which re- 
lates the triangular and hexagonal lattices as shown in 
Fig. [3J the dual lattice described in Fig. [5] is not the 
same structure of the original square lattice, since the 
location of the bond dilution/ contraction in the unit cell 
differs from the original one. Actually the extreme limit 
Pmix = 1 for the dual random planar lattice is the same 
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FIG. 8: The dual lattice of the random planar lattice depicted 
in Fig. [7] The same symbols are used as those in Fig. [JJ 
The extreme limit p m i x = 1 becomes the Union- Jack lattices, 
which is the mutually dual pair of the (4, 8 2 ) lattice. 



as the Union- Jack lattice, which is known as the mutu- 
ally dual pair of the (4,8 2 ) lattice. Our duality analysis 
can give the connection between the critical points on 
the random planar lattices with p m ix and 1 — p m ix not 
only for the extreme case. In this sense, we generalize 
the concept of the mutually dual pair into the case on 
the random planar lattice. However we cannot estimate 
the accurate locations of the critical points on this type 
of the random planar lattice only by the duality analysis, 
since the random planar lattice does not hold self-duality 
on the original square lattice as seen above. The dual- 
ity analysis can be performed only on the random planar 
lattice that satisfies the self-duality of its basic lattice as 
in Figs. d]and[Hl 



IV. ISING AND POTTS MODELS ON 
RANDOM PLANAR LATTICES 

In this section, we evaluate the critical points for the 
several spin models on the random planar lattice as in 
Fig. [2J mixture of the triangular and hexagonal lattices, 
as a benchmark test of the present method. 

Let us perform the duality analysis with the real-space 
renormalization by use of several clusters as in Fig. |H1 
By taking the cluster as the basic platform to estimate 
the location of the critical point, we can deal with the 
uniform modifications on the random planar lattice. We 
impose that Jjy on the remaining bonds without modifi- 
cation is unity in this section, in order to check that our 
method works well. In the next section, we apply our 
method to the spin glass model by introducing another 
quenched randomness on such a remaining bond Jij . 
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FIG. 9: The unit cluster for the random planar lattice s = 1 
(left panel). The right panel denotes the larger cluster s = 2 
to estimate more precise location of the critical points. The 
same symbols used as in Fig. [T] 



First we take the g-state Potts model on the random 
planar lattice. We obtain the following formula from Eq. 
([7]) by use of the four-bond cluster as in Fig. [§] (left), 
namely s — 1 cluster, to estimate the location of the 
critical points for the Potts model on the random planar 
lattice: 



2 log q = _p mix log 



q(l + (q-l)e~^ J ) 3 



1 + (g - l)e~Wc J 

+ (1 - Pmix) 

x log {(1 + (g - l)e-^ J ) 3 + (q - 1)(1 - e~^ J f] . 

(8) 

The detailed calculation is shown in Appendix [5] The 
special case with q = 2 corresponds to the critical point 
for the Ising model through T c — 2/ f3 c . In addition to 
this case, we list the critical points for the q = 3 and 
q = 4 Potts models by T c = 1/(3 C in Table Q] The duality 
relation for the mutually dual pair as confirmed in Figs. 
131 and H (p mix ,exp(-/3 c J)) ->■ (1 - p mix , (exp(-/3 c J) - 
l)/(l+(g— 1) exp(— (3 C J)), is satisfied for any value ofp m ; x 
and q. We can estimate the bond-percolation threshold 
on the random planar lattice through p* h = exp(— j3 c J) — 
1 on the random planar lattice via a technical limit of 
q — > 1 15]. The obtained results are also listed in Table 
UJ The well-known exact solutions for the triangular and 
hexagonal lattices are properly reproduced with p m ; x = 
and p m i x = 1, respectively. We confirm the Kesten's 
duality relation, which is the summation of for p mix 
and 1 — p m i x should be unity for mutually dual lattices 
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V. RANDOM-BOND ISING MODEL ON 
RANDOM PLANAR LATTICES 

Next, we apply our duality analysis to the spin glass 
model with competing interactions on the random pla- 
nar lattice given by the mixture of the triangular and 
hexagonal lattices. 

Let us perform the duality analysis on the spin glass 
model on the random planar lattice by introducing the 
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pmix 


Tc (Ising) 


T c (q = 3) T c (q = 4) 


th 
Pc 


0.00 


3.64096 


1.58492 


1.44270 


0.34730 


0.10 


3.32027 


1.44672 


1.31766 


0.37422 


0.20 


3.02041 


1.31759 


1.20097 


0.40334 


0.30 


2.74404 


1.19874 


1.09372 


0.43436 


0.40 


2.49324 


1.09104 


0.99669 


0.46680 


0.50 


2.26919 


0.99497 


0.91024 


0.50000 


0.60 


2.07194 


0.91048 


0.83425 


0.53320 


0.70 


1.90046 


0.83707 


0.76824 


0.56564 


0.80 


1.75281 


0.77389 


0.71140 


0.59666 


0.90 


1.62647 


0.71984 


0.66279 


0.62578 


1.00 


1.51865 


0.67376 


0.62134 


0.65270 



TABLE I: Critical temperatures for the Ising and several (q — 
3 and 4) Potts models, and the bond-percolation thresholds on 
the random planar lattice. The extreme limits p m i x = and 
Pmix = 1 indicate the cases on the triangular and hexagonal 
lattices, respectively. 

competing interactions on the remaining three bonds on 
the unit cluster. We restrict ourselves on the case of the 
random bond Ising model (o.e. q = 2) with the following 
Hamiltonian for simplicity 

H = -^2j ij S i S jl (9) 

where Si = ±1 is the Ising spin hereafter. The sum- 
mation is taken over all the bonds on the square lattice. 
The interaction Jy remaining after the bond dilution and 
contraction follows the bimodal distribution as 

P(J«)=^(J« + J) + (1-P)^-J)=2^ 7 , 

( 10 ) 

where p denotes the concentration of the antiferromag- 
netic interactions. Here j3 p is given by exp(— 2(3 p J) — 
p/(l — p). Notice that we can generalize our scheme to 
other distribution functions than the above bimodal in- 
teractions, such as the multi-component distribution and 
the Gaussian distribution for Jy . We evaluate the spe- 
cial critical point on the Nishimori line /3 = j3 p , the mul- 
ticritical point [ill Gil ■ We express the location of the 
multicritical point by use of the concentration of the anti- 
ferromagnetic interactions p c . The estimated values pc^ 
and p c by our duality analysis (i.e. from Eq. (J7J) are 
listed in Table HU The superscripts denote the size of the 
cluster used in evaluation of the multicritical points. The 
detailed calculation is given in Appendix. [Bj 

The value of p£ for p m i x = is consistent with the 
previous result obtained by the combination of the con- 
ventional duality and the star-triangle transformation, 
while the value obtained for p m j X = 1 is slightly different 
Q • The reason for this discrepancy on the hexagonal lat- 
tice (p m ix = 1) is that the cluster we used here is slightly 
different from that in the previous study @ . It is known 



that the results depend on the choice of the cluster in 
the case of the spin glasses As discussed in Refs. 
0, H, [IH , if we estimate the critical points by the duality 
analysis around T c of the pure Ising model, we can check 
how well the duality analysis works by comparing to the 
well-known exact solution for weak randomness . The 
values of the critical points around T c becomes closer to 
the exact solution by taking larger clusters to obtain the 
more precise estimations. 

Our duality analysis in conjunction with the gauge 
transformation and the replica method as in Ref. [20| 
can also show absence of the spin glass transition in the 
symmetric distribution (p = 1/2) on the random planar 
lattice. 

The random-bond Ising model on the random planer 
lattice has a close relationship with theoretical limita- 
tion of a class of quantum error correction codes called 
surface codes [U [2l|. Two types of the errors, namely 
bit- and phase- flip errors, on the qubits are associated 
with the antiferromagnetic interactions on the original 
and dual lattices, respectively. The frustrations of the 
system are related to the error syndromes, which are 
used to infer the locations of the errors. As naturally 
expected, if the error rate, say p, of the physical qubits is 
increased, we cannot store quantum information reliably, 
which can be viewed that increase of frustration yields 
the paramagnetic configuration. This drastic change on 
the correctability of the error in quantum error correction 
is rigorously mapped onto the phase transition from the 
ordered state to disordered one in spin glasses Thus 
the multicritical point of the random-bond Ising model 
corresponds to the theoretical limitation of the surface 
code, namely optimal error threshold, in the context of 
quantum error correction. The detailed description of 
the relation between the random-bond Ising model and 
the surface code is given in Appendix [Cl 

The results in Table |TT] for the random-bond Ising 
model correspond to the optimal error thresholds of 
the surface codes defined on the random planer lattices. 
They show fairly good agreement with those obtained 
by a suboptimal method [T^|. For instance, the sub- 
optimal method gives p c = 0.1585 for p m j X = and 
Pc = 0.0645 for p mlK — 1, while our duality analysis es- 

(2) (2) 

timated p y c — 0.16409 and p c — 0.06768, respectively. 
The discrepancies between them are around 4— 5%. No- 
tice that the suboptimal decoding can not achieve the 
optimal error thresholds which we estimate, since the 
suboptimal method estimates the critical points in the 
ground state, /3 — > oo [§, [l{|. We can also check the 
performance of the surface codes on the random planar 
lattices by computing the following quantity, 

G(p c ,p* c ) = H(p c ) + H(p* c ), (11) 

where H (p) is a binary entropy and p c and p* indicate 
the optimal thresholds for p mlx and 1 — p m i x , respectively, 
where the argument p mlx is omitted for simplicity. This 
quantity is related to the quantum Gilbert- Varshamov 
(q-GV) bound through R < 1 — G(p c ) in the limit of zero 
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asymptotic rate R -> d [H [H|. If G(p c ) does not 
exceed unity, the existence of the good error correcting 
code is guaranteed. The estimate for p m ; x = exceeds 
the q-GV bound as shown in Table [TTJ However if we use 
a larger cluster as in Fig. [5j we obtain smaller values 
for any p m i x . We expect that the further analyses by use 
of the large cluster can give more precise results, which 
would be consequently below but very close to unity as 
shown in Refs. 





pi 1 ' 




(2) 
Pc 




0.00 


0.16419 


1.00144 


0.16409 


1.00038 


0.10 


0.15323 


0.99989 


0.15308 


0.99897 


0.20 


0.14207 


0.99844 


0.14191 


0.99767 


0.30 


0.13091 


0.99727 


0.13075 


0.99656 


0.40 


0.11994 


0.99644 


0.11982 


0.99587 


0.50 


0.10940 


0.99618 


0.10930 


0.99560 


0.60 


0.09945 




0.09938 




0.70 


0.09026 




0.09018 




0.80 


0.08189 




0.08179 




0.90 


0.07437 




0.07422 




1.00 


0.06768 




0.06746 





TABLE II: Optimal error thresholds on the random planar 
lattice. The superscript expresses the cluster as in Fig. [9] 
used to estimate the threshold values. 



VI. SUMMARY 

We perform the duality analysis on the random planar 
lattices, which consist of various types of the quenched 
randomness, namely, random modification of the square 
lattice by bond dilution and contraction, and competing 
interactions on the resulting random planar lattices. The 
duality analysis with the real-space renormalization on 
the square lattice leads us to the precise locations of the 
critical points on the random planar lattice. The known 
results for the triangular and hexagonal lattices by the 
star-triangle transformation are properly reproduced in 
extreme cases of the present framework. In this sense, the 
present work is a natural extension of the star-triangle 
transformation. By applying our duality analysis, we can 
also obtain several optimal error thresholds of the surface 
codes defined on the random planar lattices. 

The realm of quantum information has been rapidly 
spread its range out. The model, which we deal with, 
is yielded from the context of quantum information (l2j 
but the method to analyze it has developed in a different 
branch of physics, statistical mechanics. We hope that 
such fascinating innovations laid across several fields of 
physics continue to be generated through various studies 
on quantum information. 
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Appendix A: Derivation of Eq. §5$ from Eq. ([?]) 

We evaluate Eq. ([7]) for the case of the g-state Potts 
model only with the single bond under the bond dilution 
and contraction on the s = 1 cluster. The renormal- 
ized Boltzmann factor in the original partition function 
is written as 




ViS(S )) 



(Al) 



where Sq is the internal spin summed over, i indicates 
each bond surrounding the internal spin Sq, n is the num- 
ber of replicas, and Vi — exp(/3Jj) — 1 depending on the 
quenched disorder. The product Yii runs over four bonds 
surrounding the internal spin. On the other hand, the 
dual renormalized Boltzmann factor is given by 



*(i) 



- k So i 



(A2) 



After equating both of the quantities as in Eq. (|7|), we 
take n — > following the recipe of the replica method. 
We reach the generic formula 



log 



(q-VYliVi + YlM + Vi) 



Zlogg. (A3) 



(9-i)+rii(i+«i) 

We set the specific bond I following the distribution 

P{Ji) = fW<5(J;) + (1 -Pmix)S(Ji - oo), (A4) 

while the remaining bonds are constant Ji = J for i =/= I. 
We then obtain 



Pmix log 



q(q + v) s 



(5-1)- 

+ (1 - Pmlx) log 
= 2\ogq. 



(1 + w) 3 
(q - l)v 3 + (q + v) 3 
(1 + v) 3 



(A5) 



This equality is reduced to Eq. ([8]) by substituting 
v — exp(/3J) — 1. If we apply the duality transforma- 
tion as v* — > q/v, we find the same formula being valid 
for p mix 1 - p mix . 
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The limit of q — >• 1 gives the bond-percolation thresh- 
olds on the random planar lattice as in the literature [l5[ . 
The leading order in e with q = 1 + e gives 



J?mix log 1 + e 



3e 



1+1! (1 + V)' c 

+ (l-Pmix)log 1 + — — + 



1+1! (l+w) 3 / 
= 21og(l + e). (A6) 

The above relation is reduced to, through p* h = v c /(l + 
v c ), 

1 - 3(1 - Pmix )pf - 3 Pmix (pf) 2 + (Pf) 3 = 0. (A7) 

This equality can be reproduced from the following exact 
solution on the general formula for the bond-percolation 
problem HH 

C(p, r, s,t) = 1 — pr — ps — rs — pt — rt — st 
+prs + prt + rst + pst = 0, (A8) 

where p, r, s and t denote the inhomogenous probabili- 
ties for the occupied bonds on the unit cell with four 
bonds of the square lattice. This equality was originally 
conjectured [25j ] from several evidence and confirmed nu- 
merically in high precisions. It was also derived in a 
different way [26| . Our case corresponds to the case with 
p = r = s = p* h and t = 1 — p m j x . This coincidence sup- 
ports to our theory for the bond-percolation thresholds 
at least. 



Appendix B: Evaluation of 

Similarly to the g-state Potts model without disorder 
in interactions, we evaluate the formula to estimate the 
location of the critical points for the case with disorder 
in interactions. The renormalized Boltzmann factor in 
the original partition function is written as 



„(!) 



*. So * * 



(Bl) 



The dual renormalized Boltzmann factor is given by 



'ils no 



- f3J 'S ) 



So i 



(B2) 



Taking n — > in Eq. ([7]), we obtain the following formula 



log 



/n 4 2cosh^J i 



2 cosh 



1 + JJtanh/3Jj 



= 2 log 2. 
(B3) 



Numerical evaluation of the above equality gives the esti- 
mation of p£ in Table [II] Similarly to the previous case, 



the bond I follows the distribution function Eq. (|A4[) . For 
other remaining bonds, the bimodal distribution func- 
tion, Eq. (|TU1) . is assigned as disorder in the interactions. 
It is straightforward to generalize the above analysis for 

(2) 

the larger cluster to obtain more precise results as p c in 
Table |TTJ Then we must take the summation over several 
internal spins and evaluate the configurational average 
for all bonds on the cluster. The complexity of the for- 
mer computation can be reduced by the special technique 
as the numerical transfer matrix II 5l. 12711 



Appendix C: Relation between spin glass model and 
quantum error correction 

Quantum information is fragile against environmental 
noise and results in the loss of coherence, that is, deco- 
herence. It is one of the main issue for the realization 
of quantum information processing to counteract the er- 
rors due to decoherence. In comparison with classical in- 
formation, quantum bit, namely qubit, is subject to not 
only bit but also phase flip errors represented by Pauli X 
and Z operators, respectively. Quantum error correction 
[28| is one of the most successful scheme to handle these 
errors, which employs multiple physical qubits to encode 
logical information into its subspace, so-called code space. 

The code space is defined by the stabilizer operators 
{Si}, which are the elements of an Abelian subgroup of 
the n-qubit Pauli group {-U,X,Y,Z}® n 0. That is, 
for any quantum state \ipL) in the code space, Si\ipL) = 
\tpL) is satisfied for all i. If an error occurs on the physical 
qubits, the state is mapped onto an orthogonal subspace 
of the code space. This means that we can detect the 
occurrence of the error by measuring the eigenvalue Si of 
the stabilizer operator Si, which we call error syndrome. 
Specifically, Si = —1 implies an existence of the error 
that anticommutes with Si. 

The surface codes Q, which have a close relationship 
with the spin glass models, are a class of quantum error 
correction codes, whose stabilizer operators are defined 
on the surfaces of lattices as follows. Suppose that a 
qubit is associated with each edge of a lattice (see Fig. 
[TOl (a)). Then, the stabilizer operators are defined for 
each face / and vertex v of the lattice as 

A f = Y[X h B V = \[Z U 

ief lev 

where C; (C = X, Z) indicates the Pauli operators of the 
qubit located on the edge I, Yiief an( ^ Yii£ V are performed 
for the qubits (edges) located around the face / and ad- 
jacent to the vertex v. The error syndromes a/ and b v 
detect Z (phase) and X (bit) errors, respectively. Note 
that if a Z error chain, say E, occurs on the code space, 
those error syndromes that are located on the boundary 
dE of the error chain E return the eigenvalue a/ = — 1 
as shown in Fig, |TD] (b). This is also the case for X errors 
and the error syndrome b v . Although we, for simplicity, 
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consider only Z errors below, the same argument works 
for X errors. 

In the error correction, we infer the location of the 
error from the error syndrome dE. Suppose that Z errors 
occur with probability pz independently for each qubit. 
Conditioned on the error syndrome dE, the probability 
of a hypothetical error chain E that has the same error 
syndrome is given by 



p(E\dE) = 



P'z 



1 



P'z 



\B\/2 




-fi(E) 



leB 



dE=dE z 



where B is the set of edges (bond) of the lattice, p' z is 
the hypothetical probability of the occurrence of Z error 
for each qubit, and 



fi(E) 



-1 for I e E 
+1 fovltiE 



which indicates the location of the error chain E. (Note 
that the actual and hypothetical error probabilities pz 
and p' z are not necessary the same.) The error correction 
is done by operating the estimated error chain E on the 
state, and hence the net effect on the code space is E+E. 
If E + E belongs to the trivial homology class as shown in 
Fig. [TU](b), the error correction succeeds (otherwise, E + 
E is a non-trivial cycle and hence acts on the code space 
non-trivially, which results in an logical error Q). By 
taking a summation over such E, we obtain the success 
probability of the error correction 



oc 



X) exp 

E\8E=dE 



where Pz is defined such that e~ fe = \Jp' z I^X — Pz)- ln 
order to simplify the summation over the same homology 
class, we introduce an Ising spin <jj 6 {+1, — 1} on each 
site i of the lattice. If E + E belongs to the trivial ho- 
mology class, there exists a configuration {cr,} such that 
fi(E) can be expressed as fi(E) — fi(E)<Ji<jj with the 
sites i and j connected by the edge I. By using this fact, 
the success probability can be reformulated as 



exp 



ft 



zZ4 



<7i<Jj 



where the edge / are replaced by the pair of sites ij , and 
fij{E) is simply denoted by jfy 

Now that the relation between the success probability 
and spin glass model becomes apparent; the success prob- 
ability of the Z error correction is nothing but the parti- 
tion function of the random-bond Ising model (RBIM), 
whose Hamiltonian is given by H = —J2{ij)Jfj a i a j- 
The location of Z errors represented by j£ = — 1 corre- 
sponds to the anti-ferromagnetic interaction due to dis- 
order, whose probability distribution is given by 

P(Jfi) = (1 - PzWfi - 1) + pzSiJ? + 1). (CI) 




frustration 
(Ising vortex) 



FIG. 10: (a) A qubit, denoted by the open circle, is located on 
each edge of the lattice. The stabilizer operators defined on 
each face and vertex, (b) Z error chains E and E, which have 
the same error syndrome. The incorrect error syndromes are 
displayed by the gray faces, (c) The associated RBIM with 
the Z error correction problem. An Ising spin is located on 
each vertex of the lattice. A coupling constant J? = fij{E) 
is assigned on each edge. The actual and hypothetical error 
chains E and E correspond to the anti-ferromagnetic interac- 
tions shown by the double lines. The excited domain wall is 
shown by the bold line. The gray-shaded area is a domain of 
Oi — — 1. The crosses indicate frustrations. Equivalently they 
are also the end points of the excited domain wall, so-called 
Ising vortexes, (d) Primal (solid lines and closed circles) and 
dual (dotted lines and open circles) lattices. According to 
the definition of the stabilizer operators shown in (a), the Z 
and X error correction problems are mapped to the RBIMs 
on the primal and dual lattices, whose concentrations of anti- 
ferromagnetic interactions are given by pz and px, respec- 
tively. 



(Recall that the hypothetical error probability p' z is given 
independently of the actual error probability pz-) The 
error syndrome af — — 1 corresponds to the frustration 
of Ising interactions. Equivalently it is also the end point 
(Ising vortex) of the excited domain wall as shown in Fig. 

nja(c). 

In order to storage quantum information reliably, the 
success probability has to be exponentially small by in- 
creasing the system size. While this is achieved with 
the error probability pz smaller than a certain value, 
so-called accuracy threshold p z , if the error probability 
is higher than it, the success probability converges to a 
constant value. This drastic change of the function p suc 
comes from a phase transition of the RBIM; in the fer- 
romagnetic phase quantum error correction succeeds. Of 
course, we should take pz = l/(e 2/3z + 1) = p' z , that 
is, the actual and hypothetical error probabilities are the 
same, in order to perform a better error correction. This 
condition corresponds to the Nishimori line [l8], on 
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the (Pz,pz) phase diagram. 

The same argument leads that the X error correction 
is mapped to the RBIM on the dual lattice (see Fig. [TO] 
(d)). The inverse temperature and concentration of the 
anti-ferromagnetic interactions are taken as flx and px 
respectively. For each surface code, we have two RBIMs 
defined on the primal and dual lattices, whose inverse 
temperatures and concentrations are Pz,x and pz.x re- 



spectively. Thus, if we modify a surface code on a primal 
lattice, both RBIMs on the primal and dual lattices are 
simultaneously changed. Specifically, the random modifi- 
cation of the stabilizer operators by a systematic removal 
of the qubits introduced in Rcf. [12( corresponds to the 
bond dilution and contraction, which have mutually dual 
actions on the primal and dual lattices. 
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